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1.  Introduction 


The  interaction  of  materials  with  magnetic  fields  has  a  long  history  in  US  Army  applications, 
such  as  magnetic  shielding  and  ferromagnetic  ally  enhanced  antennae,  but  also  more  recently  in 
the  processing  of  materials.  These  applications  require  conformable,  thin,  lightweight  materials 
with  high  magnetic  permeability,  ji.  Magnetic  permeability  is  a  macroscopic  materials  property 
defining  how  a  bulk  material  responds  to  an  external  magnetic  field  [B  =/r  H],  Effective 
magnetic  shielding  materials  are  characterized  by  1)  large  values  of  relative  magnetic 
permeability,  jir,  2)  resistance  to  local  and  global  spin-state  changes,  and  3)  isotropic  response  to 
an  external  magnetic  field.  Instead  of  blocking  an  impending  magnetic  field,  high-/rr  materials 
create  a  path  for  field  lines  to  be  deflected  around  the  shielded  area.  Anisotropy  in  //r  complicates 
the  behavior  of  shielding  materials,  as  particular  care  has  to  be  taken  to  ensure  proper  alignment 
with  an  external  field. 

Current  commercial  magnetic  shielding  materials  contain  mu-metals,  which  are  a  range  of 
nickel-iron  (Ni-Fe)  alloys  with  high  ji,  but  these  materials  are  heavy,  rigid,  and  costly. 
Furthermore,  the  current  ferromagnetic  material  being  used  in  Army  antenna  applications  is  a 
high-/!  meta-ferrite.^’^  These  materials  are  fabricated  in  limited  sizes  that  may  not  be  suitable  for 
covering  large  curved  surfaces  of  ground  or  air  vehicles.  Polymeric  matrices,  on  the  other  hand, 
are  highly  conformable,  flexible,  and  generally  light.  Polymers  have  densities  around  1-2  g/cm^, 
while  those  of  copper  (Cu)  and  mu-metal  are  almost  9  g/cm^.  Polymers  can  also  be  easily 
fabricated  to  any  size  and  shape.  By  doping  or  choosing  appropriate  monomers,  polymers  can  be 
manipulated  to  exhibit  desired  properties  such  as  mechanical  toughness  or  electrical 
conductivity. 

Recent  work  has  shown  that  adding  magnetic  particles  to  soft  materials  (carbonyl  iron  in 
polyurethane^  and  magnetite  in  skin  grafts,"^  for  example)  changes  their  bulk  response  to  a 
magnetic  field.  Doping  conductive  polyaniline  (PANI)  with  magnetite  particles  (Fe304)  creates  a 
composite  with  saturation  magnetization  approximately  100  times  that  of  the  neat  polymer.^ 
Physical  incorporation  of  magnetic  particles  in  polymer  matrices  is  well  established  and  can 
achieve  the  desired  magnetic  properties;  however,  this  approach  introduces  polymer-ceramic 
interfaces  that  can  cause  poor  dispersion  and  nonuniform  shielding.  As  an  example  of  a 
lightweight  alternative,  graphene  can  be  made  ferromagnetic  by  chemically  bonding  metal  atoms 
(both  high-  and  low-spin  ions)  onto  the  surface.^  Enhanced  ferromagnetism  in  graphene  can  be 
elicited  by  intercalating  the  sheets  with  chromium  (Cr)  ions.^  Until  now  the  emphasis  has  been 
on  Cr-based  organometallics  because  they  readily  sandwich  between  graphene  layers^  and  create 
linker  groups  between  carbon  nanotubes. ^  A  few  recent  computational  and  experimental  works 
on  organometallics  bonded  to  graphene  have  been  reported, but  the  general  trend  is  for  the 
metals’  high-spin  states  to  be  quenched  on  the  surface. 
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These  reeent  reports  on  metal-dopant-indueed  magnetism  in  graphene  are  promising 
advaneements  for  the  synthesis  of  magnetie  polymers.  Polymers  are  generally  elosed-shell, 
purely  organie  materials  (no  metal  eontent)  with  no  net  magnetie  moment.  Beeause  polymers  ean 
be  synthesized  with  aromatie  earbon  (C)  rings,  however,  similar  metal-indueed  magnetism 
should  be  possible.  Furthermore,  doping  polymers  with  magnetie  organometallie  eomplexes  may 
result  in  lightweight  low-metal-eontent  eomposites  with  good  dispersion  of  dopants  and  a 
tunable  magnetie  response. 

Prior  researeh  on  magnetic  properties  of  coordination  complexes  has  focused  on  spintronics, 
molecular  magnets,  and  quantum  computing.  In  these  applications,  it  is  desirable  to  exhibit  low 
spin-flip  energies,  i.e.,  decoupled  individual  spin  states  are  close  in  energy,  and  the  magnetic 
response  has  to  be  actively  controlled  and  maintained  via  temperature,  applied  voltage,  or  laser 
pulses.  In  contrast,  shielding  applications  require  a  passive  resistance  to  local  and  global  spin- 
state  changes.  More  specifically,  the  following  3  factors  determine  the  magnetic  permeability  of 
doped  polymer  systems: 

•  Spin  state  of  the  dopant:  Relative  magnetic  permeability,  //r,  is  related  to  spin  state,  S, 
through  the  equation 

_  ^  _  fipng^ fil(iL+SXL+S+l) 

3/cbT  ’  ^  ' 

where  juo  is  the  permeability  of  free  space,  n  is  the  concentration  of  free  spins,  g  is  the 
Lande  g  factor,  jib  is  the  Bohr  magneton,  L  is  the  orbital  angular  momentum,  ks  is  the 
Boltzmann  constant,  and  T  is  temperature.  Since  fir  is  linearly  proportional  to  S,  high  spin 
states  increase  the  interaction  with  an  external  magnetic  field. 

•  Spin-state  splitting  of  the  dopant  (or  magnetic  stabilization  energy):  This  is  the  energy 
necessary  to  change  the  spin  state.  Large  spin-state  splitting  prevents  a  decrease  in  spin 
state  under  a  (low-frequency)  magnetic  field.  A  decrease  in  spin  state  is  undesirable,  as  this 
would  subsequently  reduce  the  magnetic  character  of  the  dopant. 

•  Magnetic  coupling,  J,  between  dopants  sitting  on  different  sites  within  the  polymer  matrix: 
Positive  J  (ferromagnetic  coupling)  between  dopants  guarantees  a  constructive 
interference,  which  leads  to  high  ji.  The  strength  of  J  also  determines  the  Curie 
temperature,  above  which  spontaneous  ferromagnetic  ordering  is  lost  due  to  thermal 
fluctuations. 

This  report  describes  fundamental  research  conducted  to  assess  the  potential  for  using 
organometallie  dopants  in  polymeric  shielding  materials.  The  research  hypothesis  was  that 
ferromagnetic  polymers  can  be  realized  via  doping  with  magnetic  organometallie  complexes  and 
that  the  magnetic  structure  of  the  dopant  determines  the  bulk  magnetic  response,  i.e., 
magnetically  doped  polymeric  matrices  will  exhibit  high  magnetic  permeability  and  Curie 
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temperature  (greater  than  room  temperature).  To  test  this  hypothesis,  the  magnetie  properties  of 
several  organometallie  eomplexes  were  investigated  using  spin-polarized  quantum  ehemistry 
ealeulations.  Researeh  focused  on  organometallie  coordination  complexes  containing  high-spin 
metal  ions  (Ni,  Fe,  manganese  [Mn],  etc.)  bound  to  aromatic  ligands  with  highly  delocalized 
TT-electrons,  of  which  the  polyacenes  and  graphene  are  limiting  cases.  Such  systems  have  a  rich 
number  of  magnetic  states  based  on  metal  ion  spin  state,  ligand  chemistry,  and  the  presence  of 
counter  ions. 

This  report  describes  findings  related  to  the  aforementioned  key  factors  for  high  ji,  namely, 
preferred  spin  states  of  the  dopant  complexes  and  energy  splittings  between  different  spin  states. 
For  selected  complexes  with  high  spin  states,  the  magnetic  coupling  between  dopants  was 
evaluated  with  dimer  models,  and  the  computed  magnetic  coupling  constants,  J,  were  used  as 
input  for  classical  spin  models.  This  enabled  the  prediction  of  macroscopic  magnetic  properties 
of  the  doped  polymeric  systems  including  magnetic  susceptibility,  x  (related  to  /r),  and  Curie 
temperature,  Tc,  as  a  function  of  dopant  chemistry  and  dispersion  (defined  by  dimer  distance). 
This  report  concludes  with  recommendations  for  dopant  selection  and  general  comments  on  the 
feasibility  of  achieving  ferromagnetic  ordering  in  polymers. 


2.  Literature  Review 


The  dopants  of  interest  in  the  present  work  are  half-sandwich  organometallie  complexes 
consisting  of  a  single  transition  metal  cation  bonded  to  an  aromatic  ligand.  Each  complex 
contains  a  single  metal  atom,  M,  with  charge  +x,  and  a  ligand,  L,  with  charge  y.  Charge  neutrality 
is  imposed  by  adding  an  appropriate  number  of  CF  counter  ions  to  the  metal,  which  generates 
stool  structures  of  the  general  form  yPUC\x+y.  Such  dopants  can  be  introduced  to  polymeric 
materials  either  by  entrapment  inside  the  polymer  matrix  (physical  incorporation)  or  chemical 
incorporation  as  part  of  the  polymer  chain  backbone  or  linked  as  pendant  groups.  Many 
transition-metal  organometallie  complexes  are  paramagnetic  compounds  with  odd-electron 
counts  and  open-shell  valence  configurations. In  these  complexes,  different  spin 
configurations  result  from  the  distribution  of  unpaired  valence  electrons  among  d  orbitals  that 
have  been  split  by  the  presence  of  a  ligand  field. 

The  structures  and  magnetic  properties  of  some  transition  metal  half-sandwich  complexes 
containing  aromatic  ligands  have  been  reported  previously.  In  the  area  of  half-sandwich  and 
stool  complexes  (i.e.,  those  containing  no  more  than  one  aromatic  ligand),  the  most  extensively 
studied  ligands  are  cyclopentadienyl  [(CsHs)']'^’^’"'^  and  benzene  with  only  a  few 

works  devoted  to  cycloheptatrienyl  cyclo-octatetraene  [(CgHg)^'],^^  pyrrole 

(C4NH5),^^  and  other  heterocyclic  ligands. The  cyclopentadienyl  ligand,  in  particular,  is 
ubiquitous  in  organometallie  chemistry  and  forms  half-sandwich  complexes  with  almost  all 
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transition  metal, main-group, and  rare  earth  elements. The  strueture  and  reaetivity  of 
paramagnetic  organometallic  compounds  containing  cyclopentadienyl  are  related  to  their  spin 
state,  and  both  high-  and  low-spin  configurations  can  occur  depending  on  the  electron  count. 

Transition  metal  half-sandwich  complexes  containing  benzene  also  display  unique  magnetic 
structures.  For  example,  Pandey  et  al.^*’’^'  used  density  functional  theory  (DFT)  to  determine  the 
geometries,  energetics,  and  electronic  structures  of  neutral  and  charged  3d  transition  metal  atoms 
interacting  with  benzene.  Metal-benzene  half-sandwich  complexes  were  studied  with  and  M* 
charge  states,  and  it  was  found  that  the  spin  states  of  scandium  (Sc),  titanium  (Ti),  and  vanadium 
(V)-benzene  complexes  are  enhanced  while  those  of  Mn-,  Fe-,  cobalt-  (Co),  and  Ni-benzene 
complexes  are  reduced  from  their  values  in  the  free  metal  atom.  Only  the  magnetic  moment  of 
the  neutral  Cr  atom  remains  unchanged  when  bound  to  benzene.  These  findings  show  that  the 
spin  states  (and  resultant  magnetic  properties)  of  organometallic  complexes  can  be  tuned  by  the 
choice  of  metal-ligand  composition. 


3.  Computational  Methods 


The  organometallic  complexes  of  interest  were  half-sandwich  complexes  consisting  of  organic 
ligands  bonded  to  high-spin  first-row  transition  metals.  The  metals  considered  were  V,  Cr,  Mn, 
Fe,  Co,  and  Ni,  and  these  were  bonded  to  a  series  of  organic  ligands  containing  aromatic  C  rings, 
as  shown  in  Fig.  1.  The  ligands  were  chosen  to  probe  the  effects  of  both  ring  size  (5-8  C  atoms) 
and  ring  composition  (e.g.,  substituting  sulfur  [S]  for  C  in  thiophene  [Thio])  on  the  magnetic 
properties  of  the  final  complex. 


pyrrole  imidazole  thiophene 


Fig.  1  Organic  ligands  considered  in  this  study: 

[Cp]'  =  cyclopentadiene; 

[CFIT]^  =  cycloheptatrienyl; 

[COT]^'  =  cyclooctatetraene. 

All  organometallic  complexes  were  optimized  using  the  NWCHEM  quantum  chemistry  code,^* 
and  electronic  energies  were  calculated  with  DFT  at  the  PBE0/6-31G(d)  level  of  theory. 
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The  spin-state  splitting,  or  magnetic  stabilization  energy,  is  defined  as  the  energy  difference 
between  2  spin  states,  as  follows: 

l^^mag  =  (2) 

From  Eq.  2,  negative  stabilization  energies  indicate  that  high-spin  magnetic  ordering  is  more 
energetically  favorable  than  low-spin  ordering.  Thus  complexes  with  negative  AEmag  prefer  high 
spin  states  with  the  maximum  number  of  unpaired  valence  electrons.  Conversely,  complexes 
with  positive  AEmag  prefer  low  spin  configurations.  If  the  preferred  low  spin  state  is  exactly 
5  =  0  (i.e.,  all  valence  electrons  are  paired),  the  complex  prefers  to  be  closed  shell  and  will 
display  a  diamagnetic  response  to  an  applied  field. 

The  magnetic  coupling  constant  between  2  localized  spins  in  a  dimer  model  is  defined  by  Barone 
et  al.^^ 

j  -  ow  )  ^{.^high) 

^highi.^high'^^^ 

According  to  convention,  J>  0  indicates  ferromagnetic  exchange  while  J<  0  indicates 
antiferromagnetic  exchange.  This  convention  comes  from  the  classical  spin  Hamiltonian,  defined 
by  Heisenberg^'^  as 

H  =  -Jl.^SiSj-BY.iSu  (4) 

where  J  is  the  magnetic  coupling  between  nearest-neighbor  spins,  Si,  and  B  is  the  strength  of  the 
external  field.  With  the  defined  spin  Hamiltonian,  the  partition  function  of  a  fully  interacting 
system  of  spins  can  be  written  as  the  following: 

<2(5,r)  (5) 

from  which  all  subsequent  thermodynamic  variables  can  be  computed,  including  heat  capacity, 
magnetization,  and  magnetic  susceptibility.  Equations  for  relevant  variables  are  listed  in  Table  1. 


Table  1  Relevant  thermodynamic  variables 


Thermodynamic  Variable 

Equation 

Helmholtz  free  energy 

F{B,T)  =  -NkBT\nQ{B,T) 

Heat  capacity 

Cv  = 

fdF^ 

[dT, 

V 

Magnetization 

(  dF^ 

\dB  Jy^j, 

Magnetic  susceptibility 

/ 

z  = 

\ 

.SB  j 

T 

5 


4.  Results  and  Discussion 


4,1  Magnetic  Properties  of  Organometallic  Complexes 

Figure  2  shows  the  spin-state  splittings,  or  magnetie  stabilization  energies,  AEmag,  for 
organometallic  complexes  containing  divalent  metal  cations.  The  energies  were  calculated  with 
Eq.  2  and  are  plotted  against  atomic  number  of  the  metal  cation.  The  spin  states  and 
corresponding  DFT  energies  used  to  calculate  AEmag  are  provided  as  supplementary  data  in  the 
Appendix.  For  comparison,  stabilization  energies  calculated  with  only  Cf  ligands  are  also  plotted 
in  Fig.  2  (indicated  with  diamonds).  AFmag  for  different  cations  bonded  to  organic  ligands  spans  a 
wide  range  from  ±80  kcal/mol.  Complexes  with  large  positive  AFmag  have  a  strong  preference  for 
low  spin  states  in  which  spins  are  paired.  In  contrast,  complexes  with  negative  AFmag  prefer  high 
spin  states  with  a  majority  of  spins  unpaired. 


Impact  of  metal  ion  on  magnetic  stabilization  energy 


-•-cl 

^Cp 
■^■benzene 
-^CHT 
-t— COT 
pyrrole 
- —  imidazole 
—  thiophene 


Fig.  2  AFJmagVS.  Z  for  a  series  of  half-sandwich  organometallic  complexes  containing  divalent  metal 
cations.  AFJmag  values  calculated  with  DFT  using  the  PBEO  functional. 


Figure  3  shows  representative  results  for  3  magnetic  complexes  containing  Fe^^.  This  series 
compares  the  effects  of  ligand  composition.  Specifically,  AFmagfor  FeCpCl  is  -30  kcal/mol, 
indicating  that  this  complex  prefers  a  high  spin  state.  AFmag  becomes  stronger  (decreases  to 
-48  kcal/mol)  when  one  C  atom  in  the  5-membered  ring  is  substituted  with  nitrogen  (N)  to  form 
a  pyrrole  complex.  The  magnetic  stabilization  increases  even  more  (AFmag  decreases  to 
-56  kcal/mol)  when  a  second  C  atom  is  substituted  to  form  an  imidazole  complex.  The  increased 
magnetic  stabilization  results  from  electronic  effects  on  the  N  atoms,  and  these  results  show  that 
magnetism  in  half-sandwich  complexes  can  be  tailored  by  the  choice  of  ligand. 
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Fe[Cp]CI  Fe[Py]Cl2  Fe[lmid]Cl2 

•k 

A£'(mag)  =  -30  kcal/mol  -48  kcal/mol  -5 6  kcal/mol 

Fig.  3  Optimized  structures  and  spin  state  splittings  for  3  selected 
organometallic  complexes  containing  Permeations.  [Cp]  = 
cyclopentadiene;  [Py]  =  pyrrole;  [Imid]  =  imidazole.  A£mag 
values  calculated  with  DFT  using  the  PBEO  functional. 

The  data  in  Fig.  2  enable  a  comparison  of  the  magnetic  character  of  different  organic  ligands. 
Specially,  the  largest  spin-state  splittings  are  observed  for  inorganic  Cf  ligands,  while 
stabilization  energies  for  organic  ligands  are  comparable  across  the  range  of  atomic  numbers. 

One  notable  exception  is  the  Cr^^-benzene  complex,  which  has  a  larger  spin-state  splitting  than 
CrCb.  Aside  from  the  Cr^^-benzene  complex,  however,  the  largest  stabilization  energies  are 
achieved  with  the  Imid,  Py,  and  [COT]^'  ligands. 

Aside  from  effects  due  to  ligand,  the  choice  of  metal  cation  has  a  large  impact  on  AEmag.  When 
comparing  different  metal  ions,  there  are  some  clear  trends  in  the  data  in  Fig.  2.  The  first  is  that 
the  strongest  magnetic  stabilization  is  achieved  with  Cr^"^,  Mn^"^,  and  Fe^"^.  This  is  an  expected 
finding  given  that  these  cations  have  the  largest  number  of  unpaired  d  electrons.  In  contrast, 

Co^^,  and  display  weaker,  although  still  negative,  spin-state  splittings.  This  reveals  that  half¬ 
sandwich  complexes  containing  Co^^,  and  Ni^^  can  have  high  spin  states  but  are  less 
stabilized  with  respect  to  low  spin  states.  Thus,  high  spin-state  organometallic  complexes 
containing  Cr^"^,  Mn^^,  and  Fe^"^  are  more  favorable  for  use  in  shielding  applications  because  they 
will  display  increased  interaction  with  an  external  magnetic  field. 

Another  way  to  analyze  the  magnetic  character  of  organometallic  complexes  is  to  compare  AEmag 
with  the  magnitude  of  the  spin  state,  S.  Figure  4  shows  plots  of  AEmag  versus  2S+1  for  stable 
complexes  having  a)  different  sized  C  rings  and  b)  different  ligand  composition.  Plotting  the  data 
in  this  way  reveals  2  key  aspects  of  magnetic  structure  in  half-sandwich  complexes.  First,  high 
spin  states  (up  to  2S+1  =  5)  and  large  magnetic  stabilization  energies  are  possible  with  ligands 
containing  any  size  C  ring.  Also,  substitution  of  C  with  N  or  S  enhances  the  magnetic  properties. 
Complexes  with  Py,  Imid,  and  Thio  tend  to  be  stabilized  at  higher  spin  states  {S  =  4,5)  with 
larger  values  of  AEmag  compared  with  the  simple  5-membered  C  ring,  Cp'.  Thus,  these  ligands 
should  experience  the  largest  interaction  with  an  external  field. 
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Fig.  4  A/?mag  vs.  25+7  for  a  series  of  half-sandwich  organometallic  complexes  containing  a)  different-sized  C 
rings  and  b)  5-membered  rings  with  different  composition.  ATJmag  values  calculated  with  DFT  using  the 
PBEO  functional. 


4,2  Magnetic  Coupling  in  Dimer  Models 

The  magnetic  coupling,  J,  between  adjacent  complexes  with  high  spin  states  and  large  values  of 
AE’mag  was  determined  by  developing  dimer  models  and  computing  J  with  Eq.  3.  Figure  5  a  shows 
a  representative  dimer  model  of  2  adjacent  Fe^^-Py  complexes  separated  by  distance  d  (measured 
between  metal  centers).  Magnetic  coupling  was  calculated  at  different  dimer  distances,  and  J 
versus  d  for  dimers  of  both  FePyCb  and  FePhCb  is  plotted  in  Fig.  5b.  It  is  clear  that  magnetic 
exchange  coupling  is  positive  for  these  dimer  systems,  indicating  ferromagnetic  ordering  at 
dimer  distances  of  approximately  0. 5-1.0  nm.  However,  Jis  only  strong  at  small  dimer  distances 
and  drops  off  sharply  when  the  2  complexes  are  separated  by  more  than  1 .0  nm.  This  suggests 
that  ferromagnetic  coupling  between  complexes  on  adjacent  sites  will  be  robust  only  at  high 
dopant  levels.  If  dopants  are  dispersed  more  widely,  such  that  intercomplex  distances  are  on  the 
order  of  a  few  nanometers,  magnetic  coupling  between  complexes  will  be  quenched  unless  other 
mechanisms  for  spin  coupling  are  provided,  such  as  conjugated  polymer  backbones.  Unpaired 
electrons  on  complexes  with  high  spin  states  will  still  interact  with  an  external  magnetic  field, 
but  self-interaction  and  the  possibility  for  ferromagnetic  ordering  of  spins  will  be  negligible. 
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Fig.  5  (a)  Dimer  model  for  2  FePyCb  complexes  separated  by  distance  d.  (b)  J  vs.  d  for  dimer 
models  of  FePyCband  FePhCb.  Positive  Vindicates  ferromagnetic  coupling  between 
spins  on  adjacent  complexes. 

4,3  Macroscopic  Magnetic  Properties  from  Classical  Spin  Models 

Using  calculated  values  of  the  magnetie  exehange  interaetion  in  dimer  models,  the  maeroseopie 
magnetic  properties  of  a  3 -dimensional  (3-D)  array  of  organometallic  eomplexes  are  ealeulated 
using  the  Hamiltonian  in  Eq.  4.  The  Hamiltonian  is  solved  for  a  3-D  Ising  model  of  64,000  spins 
using  a  Monte  Carlo  FORTRAN  eode.  The  eode  was  validated  by  eomputing  magnetization,  M, 
and  heat  eapaeity,  Cv,  as  a  funetion  of  temperature  for  an  array  of  interacting  Fe  atoms  (Jpe-Fe  ~ 
16  meV)^'  under  zero  applied  field  {B  =  0).  Figure  6  shows  plots  of  M  versus  T  and  Cf  versus  T 
for  an  array  of  Fe  atoms.  As  expeeted,  the  Ising  model  prediets  ferromagnetie  ordering  and  finite 
magnetization  at  low  temperatures  (Fig.  6a).  Ordering  of  spins  due  to  magnetic  self-interaction 
weakens  at  elevated  temperatures,  as  inereasing  entropy  drives  spins  to  fiuetuate  randomly. 
Above  a  eharaeteristie  temperature — the  Curie  temperature,  Tc — spin  ordering  is  eompletely 
destroyed  and  net  magnetization  drops  to  zero  around  1,000  K.  This  is  eonsistent  with  Fe’s 
ferromagnetie-to-paramagnetic  phase  transition  at  Tc  =  1,043  K.  The  exact  value  of  Tc  predieted 
by  the  Ising  model  ean  be  estimated  from  the  plot  of  Cf  versus  T  in  Fig.  6b.  Because 
ferromagnetie-to-paramagnetie  phase  transitions  are  seeond-order  phase  transitions,  heat 
eapaeity  diverges  at  the  Curie  point.  The  predieted  value  of  Tc  from  Fig.  6b  is  850  K,  whieh  is 
fairly  elose  to  the  experimental  value  of  1,043  K.  The  value  of  Tc  for  bulk  Fe  is  underestimated 
in  this  spin  model  beeause  the  chosen  Hamiltonian  only  accounts  for  magnetic  coupling  between 
nearest  neighbors.  While  this  is  a  slight  souree  of  error  for  bulk  Fe  (magnetie  eoupling  can  occur 
between  the  seeond  and  third  nearest  neighbors),  it  is  a  reasonable  approximation  for  an  array  of 
organometallic  spins  given  that  /  drops  off  sharply  as  a  funetion  of  spin  separation  (see  Fig.  5b). 
Figure  6e  shows  a  plot  of  M  versus  T  for  the  same  array  of  Fe  atoms  but  with  inereasing 
magnitude  of  applied  magnetie  field,  B.  The  data  points  in  Fig.  6o  ean  be  used  to  ealeulate  the 
magnetie  susceptibility,  x,  by  taking  the  derivative  of  M  with  respeet  to  5  at  a  fixed  value  of  T 
(see  equation  for  /  in  Table  1). 
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Fig.  6  Plots  of  macroscopic  thermodynamic  variables  computed  for  an  interacting  Ising 
model  of  64,000  Fe  atoms  (assuming  Jpe-pe  16  meV):  a)  M  vs.  Fat  =  0,  b)  Cv 
vs.  r  at  =  0,  and  c)  M  vs.  T  for  varying  values  of  B 

The  Ising  model  was  subsequently  applied  to  a  system  of  interacting  organometallic  complexes; 
specifically,  thermodynamic  variables  were  computed  for  an  array  of  64,000  FePyCb  complexes 
separated  by  Fe-Fe  interatomic  distances  of  0.5  and  1 .0  nm.  Plots  of  relevant  variables  (at  5  =  0) 
are  shown  in  Fig.  7.  The  predicted  Curie  temperature  for  a  matrix  of  FePyCb  complexes 
separated  by  0.5  nm  is  2,100  K,  which  is  much  higher  than  that  of  bulk  Fe  because  of  the 
enhanced  magnetic  coupling  interaction,  41  meV.  In  a  real  polymer  matrix,  however,  it  is 
expected  that  dopant  loadings  will  not  be  so  dense  or  well-ordered.  As  shown  in  Fig.  5b,  the 
magnetic  coupling  drops  sharply  at  large  dimer  separations,  so  it  is  expected  that  self-interactions 
will  be  insignificant  at  sparse  loading  levels.  Therefore,  assuming  a  small  value  of  1.5  meV 
at  J  =  1.0  nm,  it  is  found  that  the  same  array  of  64,000  FePyCb  complexes  displays  no 
spontaneous  ordering  of  spins  and  effectively  zero  magnetization  at  all  temperatures. 
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T(K)  7<K) 


Fig.  7  Plots  of  macroscopic  thermodynamic  variables  computed  for  an  interacting  Ising  model  of 
64,000  FePyCh  complexes.  J  was  taken  to  be  41  meV  at  <7  =  0.5  nmand  1.5  meV  at 
(7  =  1.0  nm.  a)  M  vs.  Fat  =  0  and  b)  CfVS.  Fat  7?  =  0. 

Quantum  chemistry  calculations  show  that  magnetic  self-interactions  between  widely  dispersed 
complexes  {d>  1.0  nm)  are  quenched.  However,  complexes  with  high  spin  states  will  still 
interact  with  an  external  magnetic  field.  This  is  shown  in  Fig.  8,  which  contains  plots  of  M 
versus  T  for  organometallic  arrays  subjected  to  applied  magnetic  fields  of  varying  strength. 

As  shown  in  Fig.  8b,  the  net  magnetization  of  an  array  of  organometallic  complexes  separated  by 
1 .0  nm  is  effectively  zero  at  5  =  0  (red  line).  However,  finite  magnetization  can  be  induced  at 
low  temperatures  (i.e.,  room  temperature)  by  applying  increasingly  stronger  magnetic  fields. 
Thus,  polymers  doped  with  magnetic  organometallic  complexes  will  display  paramagnetic 
behavior  in  an  external  magnetic  field  even  if  dopants  are  widely  dispersed.  Under  these 
conditions,  the  spins  will  not  interact  with  each  other  and  will  be  “free  spins”,  responding  only  to 
thermal  fluctuations  and  external  magnetic  fields. 


T(K)  T(K) 


Fig.  8  Plots  of  M  vs.  F  for  varying  values  of  B  computed  for  an  interacting  Ising  model  of  64,000 
FePyCh  complexes  separated  by  a)  <7  =  0.5  nm  and  b)  <7  =  1.0  nm.  Red  line  corresponds  to 
fi  =  0. 
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5.  Conclusions 


Spin-polarized  quantum  chemistry  calculations  were  used  to  determine  the  structure  and 
magnetic  ordering  in  a  large  number  of  half-sandwich  organometallic  complexes.  Calculations 
validated  one  of  the  underlying  hypotheses,  which  was  that  organometallic  complexes  can  have 
high  spin  states  with  large  magnetic  stabilization  energies.  Key  chemical  design  principles  for 
achieving  high-spin  dopants  were  identified;  specifically,  the  highest  spin  states  and  largest 
magnetic  energies  are  achievable  in  dopants  with  Cr^^,  Mn^^,  and  Fe^^  bound  to  aromatic  rings. 
Moreover,  rings  substituted  with  N  and  S  (Py,  Imid,  and  Thio)  display  even  stronger  magnetic 
character  due  to  electronic  effects  on  the  ligand.  Dimer  models  were  developed  for  complexes 
having  the  largest  spin  states,  and  magnetic  coupling,  J,  was  computed  as  a  function  of  metal- 
metal  interatomic  distance.  It  was  found  that  magnetic  exchange  at  small  dimer  distances 
(~0. 5-1.0  nm)  is  large  and  positive,  indicating  ferromagnetic  ordering.  However,  J  quickly  drops 
to  zero  when  complexes  are  separated  by  more  than  1 .0  nm,  which  is  likely  to  occur  when 
dopants  are  dispersed  in  a  polymer  matrix.  Thus,  one  of  the  starting  hypotheses — namely  that 
ferromagnetic  character  could  be  induced  in  polymers — was  found  to  be  highly  dependent  on 
dopant  distance  and  not  possible  for  cases  when  complexes  are  widely  dispersed. 

A  FORTRAN  code  was  developed  to  simulate  a  3-D  classical  Ising  spin  model.  This  code  was 
validated  for  an  array  of  interacting  iron  spins  (to  mimic  bulk  Fe)  and  was  used  to  estimate  Curie 
temperature  and  magnetic  field  response  for  doped  polymeric  systems.  One  of  the  key  findings 
was  that  doped  polymeric  systems  will  behave  as  paramagnets  in  the  presence  of  an  external 
magnetic  field  even  if  magnetic  self-interactions  are  negligible,  as  would  be  the  case  when 
dopants  are  widely  dispersed.  This  indicates  that  lightweight  electromagnetic  shielding  could  be 
achievable  with  polymeric  materials  doped  with  magnetic  organometallic  complexes  even  if 
ferromagnetism  is  not  feasible.  Beyond  this  material  system,  however,  the  developed  FORTRAN 
code  is  applicable  to  any  system  of  magnetic  spins  (polymer,  metal,  ceramic,  composite,  etc.) 
and  could  be  used  to  evaluate  the  behavior  of  a  variety  of  materials  under  external  magnetic 
field. 


6.  Recommendations 


Based  on  the  results  of  this  investigation,  a  table  of  recommended  dopants  was  developed  both  to 
highlight  the  key  chemical  design  principles  and  to  guide  future  materials  selection.  Table  2  lists 
the  top  40%  of  organometallic  dopants  based  on  magnetic  performance.  All  of  these  dopants 
have  high  spin  states  with  large  values  of  AEmag  and  will  therefore  have  strong  interaction  with 
an  external  magnetic  field.  Dopants  in  bold  text  are  especially  promising  because  they  have  more 
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than  one  favorable  high  spin  state  (i.e.,  AEmag  <  0  for  multiple  values  of  S).  These  are  good 
eandidates  for  eleetromagnetie  shielding  beeause  unpaired  spins  would  align  with  the  magnetie 
component  of  the  field,  and  spin  state  flipping  could  provide  a  damping  mechanism. 

Table  2  Top  40%  of  organometallic  dopants  considered  in  the  study; 
ranking  determined  by  strength  of  magnetic  stabilization 
energy.  Bold  complexes  have  more  than  one  high  spin  state. 


Top  40% 

Top  30% 

Top  20% 

Top  10% 

Mn(Imid)Cl2 

FeHpCb 

MnCpCl 

[Fe(Thio)]2+ 

[Fe(Imid)]^^ 

MnPyCb 

CoHpCb 

IMnPhCbl^^ 

[MnCp]+ 

IMnPhCbl^- 

Cr(Imid)Cl2 

[FePy]2+ 

CrPhCh 

Mn(Thio)Cl2 

FePyCh 

[FePyCl]+ 

[Fe(Imid)Cl]+ 

Fe(Imid)Cl2 

CrHpCb 

[VCp]" 

MnHpCb 

[CrHpCb]^ 

Fe(COT) 

Co(COT) 

Mn(COT) 

MnPhCb 

Fe(Thio)Cl2 

[Fe(Thio)Cl]-" 
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Appendix.  Supplementary  Density  Functional  Theory  (DFT)  Data  for  All 

Investigated  Complexes 


This  appendix  appears  in  its  original  form,  without  editorial  change. 
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Table  A-1  Energy  difference  between  high-  and  low-spin  magnetic  states  in  transition  metal  chlorides,  M(x)Ci 


Metal  Ion 

Complex 

DFT 

Functional 

Spin  State, 

is+\ 

E  (Hartree) 

A£'mag 

(kcal/mol) 

V(II) 

VCI2 

PBEO 

1 

-1863.811784 

3 

-1863.870844 

-37.06 

B3LYP 

1 

-1864.377554 

3 

-1864.426461 

-30.69 

Cr(III) 

CrCE 

PBEO 

1 

-2424.332334 

3 

-2424.401542 

-43.43 

Cr(II) 

CrClz 

PBEO 

2 

-1964.254623 

4 

-1964.322197 

-42.40 

B3LYP 

2 

-1964.829477 

4 

-1964.883224 

-33.73 

2 

-2070.69748 

Mn(II) 

MnCE 

PBEO 

4 

-2070.76927 

-45.05 

6 

-2070.86407 

-104.54 

2 

-2643.46742 

Fe(III) 

FeCE 

PBEO 

4 

-2643.53439 

^2.02 

Fe(II) 

FeClz 

PBEO 

0 

-2183.344986 

4 

-2183.503658 

-99.57 

B3LYP 

0 

-2183.944312 

4 

-2184.086235 

-89.06 

Co(II) 

C0CI2 

PBEO 

1 

-2302.37808 

3 

-2302.51595 

-86.51 

B3LYP 

1 

-2302.981534 

3 

-2303.113368 

-82.73 

Ni(II) 

NiCE 

PBEO 

0 

-2427.92927 

2 

-2427.995758 

-41.72 

B3LYP 

0 

-2428.54387 

2 

-2428.43635 

+67.47 
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Table  A-2  Energy  difference  between  high-  and  low-spin  magnetic  states  in  transition  metal  ion-cyclopentadiene 
(Cp  )  complexes 


Metal 

Ion 

Complex 

Structure 

DFT 

Functional 

Spin 

State, 

2Y+1 

E  (Hartree) 

A£^mag 

(kcal/mol) 

V(II) 

[VCp]" 

PBEO 

1 

-1136.701991 

3 

-1136.770867 

-43.22 

B3LYP 

1 

-1137.126625 

3 

-1137.21334 

-54.41 

HSE03 

1 

-1136.847159 

3 

-1136.936618 

-56.14 

PBEO 

1 

-1597.108851 

3 

-1597.138226 

-18.43 

VCpCl 

B3LYP 

1 

-1597.723788 

3 

-1597.749817 

-16.33 

HSE03 

1 

-1597.323879 

3 

-1597.352335 

-17.86 

Cr(III) 

CrCpCb 

PBEO 

1 

-2157.61735 

3 

-2157.672298 

-34.48 

Cr(II) 

[CrCp]" 

PBEO 

2 

-1237.155921 

4 

-1237.207925 

-32.63 

B3LYP 

2 

-1237.611922 

4 

-1237.65859 

-29.28 

HSE03 

2 

-1237.325575 

4 

-1237.375202 

-31.14 

CrCpCl 

PBEO 

2 

-1697.535063 

4 

-1697.570923 

-22.50 

B3LYP 

2 

-1698.162769 

4 

-1698.190164 

-17.19 

HSE03 

2 

-1697.751124 

4 

-1697.786148 

-21.98 

Mn(II) 

[MnCp]-" 

( 

HF 

1 

-1341.808401 

A 

5 

-1341.949512 

-88.55 

C  . 

PBEO 

1 

-1343.598338 

5 

-1343.722092 

-77.66 

B3LYP 

1 

-1344.067243 

ctT 

5 

-1344.178329 

-69.71 

HSE03 

1 

-1343.770882 

5 

-1343.889657 

-74.53 

MnCpCl 

HF 

1 

-1801.650062 

c 

5 

-1801.792018 

-89.08 

PBEO 

1 

-1804.009078 

5 

-1804.094656 

-53.70 

B3LYP 

1 

-1804.64944 

VO 

5 

-1804.719071 

-43.69 

HSE03 

1 

-1804.228008 

5 

-1804.309877 

-51.37 
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Table  A-2  Energy  difference  between  high-  and  low-spin  magnetic  states  in  transition  metal  ion-cyclopentadiene 
(Cp  )  complexes  (continued) 


Metal 

Ion 

Complex 

Structure 

DFT 

Functional 

Spin 

State, 

2S+1 

E  (Hartree) 

AEmag 

(kcaFmol) 

Mn(I) 

MnCp 

PBEO 

2 

-1343.849134 

4 

-1343.892971 

-27.51 

Fe(III) 

[FeCp]2+ 

HF 

1 

- 

5 

-1454.087276 

PBEO 

1 

-1455.859648 

5 

-1455.871298 

-7.31 

B3LYP 

1 

-1456.343277 

5 

-1456.346795 

-2.21 

HSE03 

1 

-1456.036247 

5 

-1456.047116 

-6.82 

[FeCpCl]+ 

HF 

1 

-1914.015069 

5 

-1914.090913 

-47.59 

PBEO 

1 

-1916.401864 

5 

-1916.449942 

-30.17 

B3LYP 

1 

— 

5 

-1917.091622 

HSE03 

1 

-1916.658037 

FeCpCb 

5 

-1916.673819 

-9.90 

HF 

1 

— 

o  o 

5 

-2373.917495 

PBEO 

1 

-2376.788636 

5 

-2376.809397 

-13.03 

B3LYP 

1 

-2377.60434 

5 

-2377.616655 

-7.73 

HSE03 

1 

-2377.060353 

5 

-2377.080877 

-12.88 

Fe(II) 

[FeCp]-" 

c 

HF 

0 

-1454.445351 

4 

-1454.527657 

-51.65 

c 

PBEO 

0 

-1456.318156 

4 

-1456.371941 

-33.75 

B3LYP 

0 

-1456.800703 

FeCpCl 

cT 

4 

-1456.844734 

-27.63 

HSE03 

0 

-1456.491531 

4 

-1456.542533 

-32.00 

HF 

0 

-1914.284108 

c 

4 

-1914.383227 

-62.20 

PBEO 

0 

-1916.690788 

4 

-1916.749366 

-36.76 

cipyc 

B3LYP 

0 

-1917.342553 

4 

-1917.389757 

-29.62 

HSE03 

0 

-1916.911793 

Vi/ 

4 

-1916.967178 

-34.75 

22 


Table  A-2  Energy  difference  between  high-  and  low-spin  magnetic  states  in  transition  metal  ion-cyclopentadiene 
(Cp  )  complexes  (continued) 


Metal 

Ion 

Complex 

Structure 

DFT 

Functional 

Spin 

State, 

2S+1 

E  (Hartree) 

AEmag 

(kcaFmol) 

PBEO 

1 

-1575.364936 

3 

-1575.357752 

B3LYP 

1 

-1575.85711 

4.51 

[CoCp]+ 

3 

-1575.845496 

HSE03 

1 

-1575.544735 

7.29 

Co(II) 

3 

-1575.530186 

PBEO 

1 

-2035.734083 

9.13 

3 

-2035.738356 

CoCpCl 

B3LYP 

1 

-2036.393756 

-2.68 

3 

-2036.393183 

HSE03 

1 

- 

0.36 

3 

-2035.957733 

Ni(III) 

NiCpC12 

1 

-2621.315467 

PBEO 

3 

-2621.290785 

15.49 

HF 

0 

-1698.895451 

c 

2 

-1698.936136 

-25.53 

PBEO 

0 

-1700.856128 

[NiCp]+ 

2 

-1700.896852 

-25.55 

B3LYP 

0 

-1701.361268 

2 

-1701.396748 

-22.26 

HSE03 

0 

-1701.031379 

Ni(II) 

2 

-1701.070874 

-24.78 

HF 

0 

-2158.745278 

c 

2 

-2158.786974 

-26.16 

PBEO 

0 

-2161.222378 

NiCpCl 

2 

-2161.269587 

-29.62 

B3LYP 

0 

-2161.893686 

2 

-2161.936897 

-27.12 

HSE03 

0 

-2161.444688 

2 

-2161.490578 

-28.80 
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Table  A-3  Energy  difference  between  high-  and  low-spin  magnetic  states  in  transition  metal  ion- 
benzene  (Ph)  complexes 


Metal  Ion 

Complex 

DFT 

Functional 

Spin 

State, 

2S+\ 

E  (Hartree) 

A£^mag 

(kcal/mol) 

V(II) 

[VPhCF]- 

PBEO 

2 

4 

-2556.01914 

-2556.04215 

-14.44 

2 

-2095.87438 

V(II) 

VPhCE 

PBEO 

4 

-2095.88844 

-8.82 

6 

-2095.77056 

65.15 

V(IV) 

[VPhCb]-" 

PBEO 

2 

4 

-2555.67583 

-2555.6336 

26.50 

V(V) 

[VPhCE]^^ 

PBEO 

1 

3 

-2555.1926 

-2555.16097 

19.85 

Cr(II) 

[CrPhCF]- 

PBEO 

1 

3 

-2656.39646 

-2656.44955 

-33.31 

2 

-2196.20995 

Cr(II) 

CrPhCE 

PBEO 

4 

-2196.30506 

-59.68 

6 

-2196.33057 

-75.69 

2 

-2656.35458 

Cr(III) 

CrPhCF 

PBEO 

4 

-2656.39349 

-24.42 

6 

-2656.33357 

13.18 

Cr(V) 

[CrPhCF]^^ 

PBEO 

1 

3 

-2655.55977 

-2655.55064 

5.73 

1 

-2762.81288 

Mn(I) 

[MnPhCb]^- 

PBEO 

3 

-2762.85041 

-23.55 

5 

-2762.91841 

-66.22 

1 

-1842.67306 

Mn(I) 

MnPhCl 

PBEO 

3 

-1842.68342 

-6.50 

5 

-1842.73151 

-36.68 

2 

-2302.76584 

Mn(II) 

MnPhClz 

PBEO 

4 

-2302.79628 

-19.10 

6 

-2302.85177 

-53.92 

2 

-2762.49734 

Mn(IV) 

[MoPhCb]"" 

PBEO 

4 

-2762.55561 

-36.57 

6 

-2762.53437 

-23.24 

1 

-2761.95828 

Mn(V) 

[MnPhCb]^^ 

PBEO 

3 

-2762.03255 

-46.60 

5 

-2762.03395 

-47.48 

1 

-2415.44017 

Fe(II) 

FePhCb 

PBEO 

3 

-2415.46285 

-14.23 

5 

-2415.50193 

-38.76 

Fe(III) 

FePhCb 

PBEO 

2 

4 

-2875.46512 

-2875.50714 

-26.36 

2 

-2534.4861 

Co(II) 

CoPhCb 

PBEO 

4 

-2534.50532 

-12.06 

6 

-2534.39226 

58.88 
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Table  A-3  Energy  difference  between  high-  and  low-spin  magnetic  states  in  transition  metal  ion- 
benzene  (Ph)  complexes  (continued) 


Metal  Ion 

Complex 

DFT 

Functional 

Spin 

State, 

2S+\ 

E  (Hartree) 

A£^mag 

(kcal/mol) 

1 

-2659.97365 

Ni(II) 

NiPhCE 

PBEO 

3 

-2659.99398 

-12.76 

5 

-2659.87838 

59.78 

2 

-3119.99592 

Ni(III) 

NiPhCE 

PBEO 

4 

-3120.01725 

-13.39 

6 

-3119.88663 

68.58 
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Table  A-4  Energy  difference  between  high-  and  low-spin  magnetic  states  in  transition  metal  ion- 
cycloheptatrienyl  (Hp"^)  complexes 


Metal  Ion 

Complex 

DFT 

Functional 

Spin 

State, 

is+\ 

E  (Hartree) 

A£^mag 

(kcal/mol) 

2 

-2594.57251 

V(II) 

VHpCF 

PBEO 

4 

-2594.51416 

36.61 

1 

-2694.90808 

Cr(II) 

CrHpCF 

PBEO 

3 

-2694.98643 

-49.16 

5 

-2694.99086 

-51.95 

2 

-2694.65958 

Cr(III) 

[CrHpCb]^ 

PBEO 

4 

6 

-2694.749 

-56.11 

2 

-2801.40881 

Mn(II) 

MnHpCb 

PBEO 

4 

-2801.48713 

-49.15 

6 

-2801.52384 

-72.18 

1 

-2914.06604 

Fe(II) 

FeHpCb 

PBEO 

3 

-2914.07081 

-2.99 

5 

-2914.17492 

-68.32 

2 

-2913.77955 

Fe(III) 

[FeHpCb]^ 

PBEO 

4 

-2913.76618 

8.39 

2 

-3033.08797 

Co(II) 

CoHpCb 

PBEO 

4 

-3033.17295 

-53.32 

6 

-3033.10229 

-8.99 

1 

-3158.62106 

Ni(II) 

NiHpCb 

PBEO 

3 

-3158.61129 

6.13 

5 

-3158.59114 

18.77 

2 

-3158.35824 

Ni(III) 

[NiHpCb]^ 

PBEO 

4 

6 

-3158.26236 

60.16 
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Table  A-5  Energy  difference  between  high-  and  low-spin  magnetic  states  in  transition  metal  ion-cyclooctatetraene 
([COT]^  )  complexes 


Metal  Ion 

Complex 

Structure 

DFT 

Functional 

Spin 

State, 

2Y+1 

E  (Hartree) 

\Emag 

(kcal/mol) 

V(II) 

V(COT) 

PBEO 

1 

-1252.90288 

3 

-1252.949255 

-29.10 

B3LYP 

1 

-1253.493913 

3 

-1253.532379 

-24.14 

HSE03 

1 

-1253.151167 

3 

-1253.1771 

-16.27 

Cr(III) 

Cr(COT)Cl 

PBEO 

1 

-1813.470457 

3 

-1813.503411 

-20.68 

Cr(II) 

Cr(COT) 

PBEO 

2 

-1353.340621 

4 

-1353.385175 

-27.96 

B3LYP 

2 

-1353.939721 

4 

-1353.977031 

-23.41 

HSE03 

2 

-1353.572457 

4 

-1353.615258 

-26.86 

Mn(II) 

Mn(COT) 

PBEO 

1 

-1459.766615 

5 

-1459.88146 

-72.07 

B3LYP 

1 

-1460.412796 

5 

-1460.481068 

-42.84 

HSE03 

1 

-1460.034791 

5 

-1460.113929 

-49.66 

Fe(III) 

[Fe(COT)]+ 

HF 

1 

- 

5 

-1569.914411 

PBEO 

1 

-1572.296859 

5 

-1572.288777 

5.07 

Cj 

B3LYP 

1 

-1572.920006 

r 

5 

-1572.922376 

-1.49 

A 

HSE03 

1 

-1572.5329 

5 

-1572.540011 

-4.46 

Fe(COT)Cl 

PBEO 

1 

-2032.616565 

5 

-2032.663545 

-29.48 

B3LYP 

1 

-2033.429486 

5 

-2033.446452 

-10.65 

HSE03 

1 

-2032.926909 

5 

-2032.948066 

-13.28 
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Table  A-5  Energy  difference  between  high-  and  low-spin  magnetic  states  in  transition  metal  ion-cyclooctatetraene 
([COT]^  )  complexes  (continued) 


Metal  Ion 

Complex 

Structure 

DFT 

Functional 

Spin 

State, 

2Y+1 

E  (Hartree) 

\Emag 

(kcal/mol) 

Fe(II) 

Fe(COT) 

PBEO 

0 

-1572.458193 

4 

-1572.540577 

-51.70 

B3LYP 

0 

-1573.085564 

4 

-1573.155745 

-44.04 

HSE03 

0 

-1572.69417 

4 

-1572.773865 

-50.01 

Co(II) 

Co(COT) 

PBEO 

1 

-1691.480649 

3 

-1691.54927 

-43.06 

B3LYP 

1 

-1692.163048 

3 

-1692.179126 

-10.09 

HSE03 

1 

-1691.766314 

-1691.784899 

-11.66 

Ni(III) 

Ni(COT)Cl 

PBEO 

1 

-2277.177861 

3 

-2277.171639 

3.90 

Ni(II) 

Ni(COT) 

PBEO 

0 

-1816.923405 

2 

-1816.803252 

75.40 

B3LYP 

0 

-1817.568768 

2 

-1817.696978 

-80.45 

c 

HSE03 

0 

-1817.307678 

2 

-1817.289441 

11.44 
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Table  A-6  Energy  difference  between  high-  and  low-spin  magnetic  states  in  transition  metal  ion-pyrrole  (Py) 
complexes 


Metal 

Ion 

Complex 

Structure 

DFT 

Functional 

Spin 

State, 

2S+1 

£■  (Hartree) 

A£'mag 

(kcal/mol) 

Cr(II) 

CrPyCb 

PBEO 

2 

-2174.256742 

4 

-2174.293651 

-23.16 

Mn(II) 

MnPyCE 

PBEO 

1 

-2280.715037 

5 

-2280.82386 

-68.29 

Fe(II) 

[FePy]2+ 

c 

HF 

0 

-1470.644221 

4 

-1470.740657 

-60.51 

PBEO 

0 

-1472.500893 

MB 

4 

-1472.575717 

-46.95 

B3LYP 

0 

-1472.997941 

4 

-1473.062156 

-40.30 

HSE03 

0 

-1472.677264 

0 

4 

-1472.741404 

-40.25 

[FePyCl]+ 

HF 

0 

A 

-1930.652465 

PBEO 

0 

-1933.056694 

C7X 

4 

-1933.12817 

-44.85 

B3LYP 

0 

-1933.724937 

4 

-1933.783817 

-36.95 

HSE03 

0 

-1933.284878 

4 

-1933.35085 

-41.40 

FePyCb 

PBEO 

0 

-2393.39154 

2 

-2393.421948 

-19.08 

xJWCn 

4 

-2393.474424 

-52.01 

r  Ci 

B3LYP 

0 

-2394.221505 

4 

-2394.297869 

-47.92 

J 

HSE03 

0 

-2393.66438 

4 

-2393.745053 

-50.62 
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Table  A-7  Energy  difference  between  high-  and  low-spin  magnetic  states  in  transition  metal  ion-imidazole  (Imid) 
complexes 


Metal 

Ion 

Complex 

Structure 

DFT 

Functional 

Spin 

State, 

2S+\ 

£■  (Hartree) 

A£'mag 

(kcal/mol) 

Cr(II) 

Cr(Imid)Cl2 

PBEO 

2 

-2190.273369 

4 

-2190.356811 

-52.36 

Mn(II) 

Mn(Imid)Cl2 

PBEO 

1 

-2296.748122 

5 

-2296.879117 

-82.20 

Fe(II) 

[Fe(Imid)]2+ 

PBEO 

0 

-1488.491421 

4 

-1488.618035 

-79.45 

B3LYP 

0 

-1489.000355 

4 

— 

HSE03 

0 

-1488.670679 

4 

— 

[Fe(Imid)Cl]-" 

PBEO 

0 

-1949.06401 

c 

4 

-1949.179728 

-72.61 

> 

B3LYP 

0 

-1949.741627 

Jt 

4 

-1949.850607 

-68.39 

HSE03 

0 

-1949.292257 

Fe(Imid)Cl2 

4 

-1949.40474 

-70.58 

PBEO 

0 

-2409.406799 

4 

-2409.502944 

-60.33 

ctT 

B3LYP 

0 

-2410.249549 

4 

-2410.338716 

-55.95 

HSE03 

0 

-2409.682977 

4 

-2409.776905 

-58.94 
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Table  A-8  Energy  difference  between  high-  and  low-spin  magnetic  states  in  transition  metal  ion- 
thiophene  (Thio)  complexes 


Metal  Ion 

Complex 

DFT 

Functional 

Spin 

State, 

is+\ 

E  (Hartree) 

A£^mag 

(kcal/mol) 

Cr(II) 

Cr(Thio)Cl2 

PBEO 

2 

-2516.994806 

4 

-2517.028793 

-21.33 

Mn(II) 

Mn(Thio)Cl2 

PBEO 

1 

-2623.455421 

5 

-2623.554667 

-62.28 

Fe(II) 

[Fe(Thio)]2+ 

PBEO 

0 

-1815.226675 

4 

-1815.302382 

^7.51 

[Fe(Thio)Cl]+ 

PBEO 

0 

-2275.792916 

4 

-2275.858054 

-40.87 

Fe(Thio)Cl2 

PBEO 

0 

-2736.129101 

4 

-2736.208693 

-49.94 
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List  of  Symbols,  Abbreviations,  and  Acronyms 


c 

carbon 

Co 

cobalt 

COT 

cyclo-octatetraene,  (CgHg) 

Cp 

cyclopentadienyl,  (CsHs)' 

Cr 

chromium 

Cu 

copper 

DFT 

density  funetional  theory 

Fe 

iron 

Hp 

oyeloheptatrienyl,  (CtH?)^ 

Imid 

imidazole,  C3N2H4 

Mn 

manganese 

N 

nitrogen 

Ni 

niekel 

PANI 

polyaniline 

Ph 

benzene,  CeHe 

Py 

pyrrole,  C4NH5 

S 

sulfur 

Sc 

seandium 

Thio 

thiophene,  C4SH4 

Ti 

titanium 

V 

vanadium 

Nomenclature 

Z 

atomie  number 

s 

spin  state 
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AEmag  magnetic  stabilization  energy 

J  magnetic  coupling  constant 

T  temperature 

Cv  heat  capacity 

B  magnetic  field 

M  magnetization 

X  magnetic  susceptibility 

fi  magnetic  permeability 

jir  relative  magnetic  permeability 

jiQ  permeability  of  free  space,  47r  x  1 0'^  N/A^ 

jiB  Bohr  magneton,  9.274  x  1 0'^"^  J/T 
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